Mon. Not. R. Astron. Soc. 000, 1-11 (0000) Printed 2 February 2008 (MN WT^X style file vl.4) 



Spatially resolved stellar populations in the isolated 
elliptical NGC 821 

Robert N. Proctor^ Duncan A. Forbes^, Amy Forestell^ & Karl Gebhardt^ 

^ Centre for Astrophysics & Supercomputing, Swinburne University, Hawthorn VIC 3122, Australia 
Email: rproctor@astro.swin.edu.au, dforbes@astro.swin.edu.au 

^ Department of Astronomy, University of Texas at Austin, CI4OO, Austin, TX 78712 
Email: amydove@astro.as.utexas.edu, gebhardt@astro. as.utexas.edu 



2 February 2008 



ABSTRACT 

We present the analysis of Lick absorption-line indices from three separate long-slit 
spectroscopic observations of the nearby isolated elliptical galaxy NGC 821. The three 
data sets present a consistent picture of the stellar population within one effective 
radius, in which strong gradients are evident in both luminosity-weighted age and 
metallicity. The central population exhibits a young age of ~ 4 Gyr and a metallicity ~ 
3 times solar. At one effective radius the age has risen to ~ 12 Gyr and the metallicity 
fallen to less than ~ i solar. The low metallicity population around one effective radius 
appears to have an exclusively red horizontal branch, with no significant contribution 
from the blue horizontal branch evident in some globular clusters of the same age 
and metallicity. Despite the strong central age gradient, we demonstrate that only a 
small fraction (<10%) of the galaxy's stellar mass can have been created in recent star 
formation events. We consider possible star formation histories for NGC 821 and find 
that the most likely cause of the young central population was a minor merger or tidal 
interaction that caused NGC 821 to consume its own gas in a centrally concentrated 
burst of star formation 1 to 4 Gyr ago. 
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1 INTRODUCTION 

The advent of the Lick system of spectral indices (Worthey 
1994) has provided researchers with a useful tool in the dis- 
entanglement of age and metallicity effects in the integrated 
spectra of stellar populations. This tool has already been 
used to uncover interesting trends in the properties of the 
stellar populations in the centres of galaxies (e.g. Worthey, 
Faber & Gonzalez 1992; Gorgas et al. 1997; Greggio 1997; 
Trager et al. 2000; Proctor & Sansom 2002; Terlevich & 
Forbes 2002; Proctor et al. 2004b). However, few studies to 
date have probed significantly beyond the central regions. 
Most studies have hence only sampled a small, and perhaps 
biased, fraction of the galaxy mass. The few studies that 
probe off-centre regions (e.g. Mehlert et al. 2003; Proctor 
2002, Sanchez-Blazquez et al. 2005) generally find metallic- 
ity gradients to be strong but both age and [Mg/Fe] gra- 
dients to be small or non-existent. Here we detail the ap- 
plication of Lick indices to the spatially resolved integrated 
spectra of NGC 821, which we show to possess strong gra- 
dients in both metallicity and age. 

NGC 821 is a rare example of an isolated elliptical 
galaxy. At m-M = 31.91 (Tonry et al. 2001) it is the closest 



early-type galaxy in the isolated galaxy sample of Reda et 
al. (2004). Reda et al. showed that NGC 821 conforms to 
the known scaling relations for elliptical galaxies (i.e. the 
Fundamental Plane and colour- magnitude relation). It has 
a moderate luminosity {Mv ~ -21.12) and highly elongated 
isophotes (E6). An isophotal analysis indicates the presence 
of an edge-on stellar disk along the major axis (e.g. Goud- 
frooij et al. 1994; Scorza & Bender 1995; de Souza, Gadotti 
& dos Anjos 2004). It also reveals disk-like kinematics (Em- 
sellem et al. 2004). Otherwise the galaxy shows no morpho- 
logical fine structure (Michard & Prugniel 2004; Goudfrooij 
et al. 1994). In a recent catalogue of galaxy ISM proper- 
ties, Bettoni, Galletta & Garca-Burillo (2003) measured a 
dust mass that is at the very low end of the distribution 
for elliptical galaxies (Forbes 1991), while HI and molecular 
hydrogen are undetected. The ISM properties of NGC 821 
are therefore consistent with stellar mass loss alone and do 
not require any external origin. 

In previous studies, the luminosity-weighted age of the 
stellar population in the central region of NGC 821 has 
been measured to be 4.2 Gyr (Denicolo et al. 2005), 7.7 Gyr 
(Trager et al. 2000), 11.5 Gyr (Caldwell, Rose & Concannon 
2003) and 12.5 Gyr (Vazdekis, Trujillo & Yamada 2004). 
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Thus considerable variation exists in the hterature from 
claims of very young to very old ages. 

In this work wc use Lick indices to determine the ra- 
dial profiles in age, metallicity and 'a'-element abundance 
ratios in NGC 821 out to one effective radius (Ire; the ra^ 
dius within which half the galaxy light is contained). We use 
the x^-minimisation technique of comparing a large number 
of indices to SSP models to derive ages and motallicitics 
as outlined in Proctor & Sansom (2002). We find a natural 
explanation for the range of ages measured in previous stud- 
ios. Wc also briefly discuss possible evolutionary histories for 
NGC 821. 



2 DATA 

In this paper long-slit spectroscopy of three separate studies 
of NGC 821 are analysed. Each data set covered a wave- 
length range suitable for the measurement of Lick indices. 
Trager et al. (1998) index definitions are used throughout. 

The first data set was obtained using GMOS on Gem- 
ini (North). These data were collected using a 0.5 arcsec 
slit positioned along the major axis of the galaxy. The ob- 
servations covered the wavelength range 4100 - 6930 A at 
3.5 A resolution and a dispersion of 1.77 A. Three 30 minute 
exposures were obtained. The data were reduced using the 
Gemini pipe-line. Residuals to the wavelength calibration 
solution were <0.1 A. The three exposures were co-added 
prior to the extraction of spectra. During extraction, the 
spectra were spatially binned to provide a signal-to-noise 
>20 at all radial positions. Spectra from the central and 
outer regions are shown in Fig. 1. Reliable index estimates 
were obtained for 18 Lick indices from Ca4227 to Fe5782. 
The GMOS observations were obtained as part of the au- 
thor's program to measure indices out to large radii (~2re). 
Unfortunately, only 1.5 of the 8 hours requested exposure 
time was actually observed. This constrains the index anal- 
ysis to within ~ Ire. Neither Lick-index nor flux-calibration 
standard stars were observed. The spectra could therefore 
not be flux calibrated by the normal means, and the indices 
could not be fully calibrated to the Lick system. The effects 
of these uncertainties will be discussed in Section 2.1. 

The second data set was obtained using the Hobby- 
Ebcrly Telescope (HET). These data were originally col- 
lected for the purposes of analysis of the spatially resolved 
stellar kinematics (Forestall & Gebhardt 2005). The data 
include long-slit spectroscopy along both major and minor 
axes using a 1 arcsec slit. The observations covered the wave- 
length range 4300 - 7250 A at 4.5 A resolution and a disper- 
sion of 1.95 A. This wavelength range allowed the accurate 
measurement of 15 Lick indices from Fe4383 to TiOi (NaD 
was excluded due to the effects of interstellar absorption). 
The data were reduced and binned as described in Forestall 
& Gebhardt (2005), but were not fiux calibrated. A signal- 
to-noise suitable for index determination was achieved out 
to approximately Ire. As with the Gemini data, Lick index 
standard stars were not observed so the indices could not be 
fully calibrated to the Lick system (see Section 2.1). 

The final data set was obtained from the study of 
Sanchez-Blazquez et al. (2005; hereafter SB05) as line in- 
dex measurements. The spectrum from the central region 
was also supplied for the purposes of calibration of the other 
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Figure 1. Gemini spectra from the central (upper) and outer 
(lower) regions of NGC 821 are shown. The spectra have signal- 
to-noise of 80 and 20 respectively. The presented spectra have 
been broadened to the Lick resolution of the H/3 index (8.5 A) by 
convolution with a Gaussian. An interpolation across a chip gap 
is evident at ~5000 A. 

two data sets (see Section 2.1). The data were collected us- 
ing the 3.5 m telescope at the German-Spanish Astronomical 
Observatory at Calar Alto. Long-slit spectroscopy was car- 
ried out with the 2.1 arcsec slit aligned to the parallactic 
angle, i.e. along a line at 25° to the major axis. The obser- 
vations covered the wavelength range 3570 - 5770 A at 2.6 A 
resolution and a dispersion of 1.08 A. This allowed accurate 
measurement of 19 Lick indices from H(5^ to Fe5335. Al- 
though these data were of significantly lower signal-to-noise 
than the Gemini and HET data, the spectra were accurately 
flux calibrated prior to index measurements and fully cali- 
brated to the Lick system. 



2.1 Calibrations 

The calibration of Lick indices is a two stage process. The 
first is to compensate for the effects of internal velocity dis- 
persion broadening. For the Gemini data velocity disper- 
sions were estimated by cross-correlation of the galaxy spec- 
tra with stellar spectra observed with the same instrument 
configuration. The IRAF command fxcor was used. For the 
HET data set measurements of kinematics were taken from 
Forestall & Gebhardt (2005). For both Gemini and HET 
data sets indices were measured and velocity dispersion cor- 
rected to the appropriate resolution as described in Proctor 
& Sansom (2002). For the SB05 data set, both velocity- 
dispersion corrected indices and kinematics were supplied by 
the original authors. A comparison of the radially resolved 
velocity dispersions from the three studies is shown in Fig. 2. 
Data points are shown at their projected distance from the 
galaxy centre in this and all other radial plots in this work. 
Despite the differing slit widths and orientations between 
the three studies, the comparison in Fig. 2 can be seen to 
be generally good. Agreement is also good with the central 
velocity dispersions quoted in the literature (e.g. Bender, 
Sagha & Gerhard 1994; 208 km s"\ Di NeUa et al. 1995; 
212 km s"\ Denicolo et al. 2004; 212 km s"^). 

The second stage of the calibration to the Lick sys- 
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Figure 2. Radial profiles of velocity dispersion from all three 
data sets. An effective radius of 50 arcsec (dc Vaucouleurs ct al. 
1991) is assumed tlirougliout tliis work. Gemini data are shown 
as filled stars, the HET data as filled circles and the SB05 data as 
filled squares. The agreement between studies is generally good. 



tern is to measure offsets between the indices observed in 
a number of Lick library stars and their published values. 
This process primarily compensates for differences in flux- 
calibration between the observations and the stars used to 
construct the SSP models. The SB05 data sot was fully cal- 
ibrated in this way (using 20 standard stars) and the aver- 
age offsets supplied with the indices. As noted above, the 
SB05 spectra were also accurately flux-calibrated. The cen- 
tral spectrum supplied with this data set was therefore used 
to estimate flux-calibration curves for both Gemini and HET 
data sets. This was achieved by third-order continuum fit- 
ting of the ratio of the fluxed central spectrum from SB05 
and the un-fluxed central Gemini and HET spectra. The 
flux-calibration curves generated were then applied to their 
respective data sots at all radii. Indices wore measured in 
both flux-calibrated and un-flux-calibrated spectra for the 
purposes of comparison (see Section 3). We shall refer to 
these as the fluxed and un-fluxed indices in the Gemini and 
HET data sets. However, it should be noted the the SB05 
spectra did not cover the wavelength range of indices red- 
der than Fe5335. The un-fluxed values of these indices were 
therefore used throughout this work. Average errors and the 
index offsets caused by the application of the flux calibra- 
tion curve are given in Table 1. Both flux and Lick offsets 
are generally smaller than the average index errors. 

Since our analysis probes the low surface brightness 
outer regions of the galaxy we must consider the effects of 
scattered light on our index determinations. For the Gemini 
data this can be achieved by consideration of count rates 
in regions of the slit that are masked off (there is one re- 
gion at each end of the slit and two approximately one-third 
and two-thirds along the length of the slit). Modelling these 
count rates as scattered light with a Gaussian distribution, 
it is possible to show that no more than 10% of the light 
in the most affected regions (the outermost points) can be 
from scattered light. Such contamination levels would re- 
sult in the systematic underestimation of the indices by an 
amount of up to 10% in the outer regions. However, since 
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Figure 3. Radial profiles of twelve Lick indices are shown. Filled 
symbols are as Fig. 2. The central values of Denicolo et al. (2004) 
(open circles) and Gonzalez (1993) (open squares) are also shown. 
Agreement between studies is generally good, with the exception 
of Mgi and Mg2 (see text). Gradients are evident in all metallicity 
sensitive indices, while H/3 shows no gradient. 



this is an upper limit on the contamination levels no attempt 
has been made to correct indices for this affect. The impact 
of this uncertainty on our derived ages and metallicities will 
be discussed in Section 4. 



2.2 Radial profiles in indices 

Radial profiles of twelve indices (which have been flux- 
calibrated but not corrected to the Lick system) are shown 
in Fig. 3. Agreement between the three studies is gener- 
ally good, despite the differing slit widths and orientations. 
Mgi and Mg2 arc notable exceptions to this good agreement 
(note the data shown are the values from the flux calibrated 
spectra). The cause of this discrepancy is not known, but is 
probably associated with poor flux-calibration of the Gemini 
and HET data sets. However, this problem does not effect 
our overall conclusions as the method employed permits the 
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Figure 4. Plot of H/3 against [MgFe]. Open circles are the Galac- 
tic globular clusters of Proctor, Forbes & Beasley (2004a). Open 
squares are extra-galactic sub-populations from the meta-study 
of Strader et al. (2005). Remaining symbols are as Fig. 2. Grid 
lines show the SSP models of TMB03 (solid lines) and TMK04 
(dashed lines) for an 'a'— element enhancement of -|-0.3 dex. The 
near horizontal lines show ages; from bottom to top of 15, 12, 10, 
8, 5, 3, 2, 1 Gyr. The thick lines denote the oldest age grid lines. 
The near vertical lines show metallieity; from right to left of -1-0.4, 
-1-0.25, 0.0 -0.25, -0.5, -0.75, -1.0, -1.25, -1.5, -1.75, -2.0 and -2.25 
dex. TMB03 models are those for populations with an exclusively 
red horizontal branch. The poor agreement of the TMK04 SSP 
models (which include a blue horizontal branch component) with 
the low metallieity population in NGC 821 is clearly evident. Both 
sets of models suggest that radial gradients are present in both 
age and metallieity. 

exclusion of such troublesome indices (see Section 3) . Nega- 
tive gradients are evident in all metallieity sensitive indices. 
The H/3 index, on the other hand, shows no significant grar 
dient. 

Instructive index-index plots are shown in Figs 4 and 
5*. Fig. 4 shows that in comparison to both Thomas, Maras- 
ton & Bender (2003; hereafter TMB03) and Thomas, Maras- 
ton & Korn (2004; hereafter TMK04) models, the galaxy ap- 
pears to possess both age and metallieity gradients. This fig- 
ure also demonstrates that the central indices in this galaxy 
exceed those modelled by TMB03 or TMK04. It was there- 
fore necessary to linearly extrapolate the models to an over- 
all metallieity ([Z/H]) of +0.8 dex in order to fit these data. 
It is interesting to note that the low metallieity data points 
(with [MgFe]<2.0 A) show much better agreement with the 
red horizontal branch models of TMB03 (solid lines in Fig. 
4) than the models of TMK04 (dashed lines) which include a 
blue horizontal branch component. Indeed, the galaxy data 
points appear to follow the locus of those Galactic globular 
clusters which possess exclusively red horizontal branch mor- 
phologies (Proctor et al. 2004a). Fig. 5 also shows the consis- 
tency of the data with the SSP models of TMB03 with 'a'- 
element enhancements of -1-0.3 dex, indicating that NGC 821 
possesses a similar abundance ratio to other bright elliptical 

* [MgFe] = ^Mgb x (Fe5270 + Fe5335)/2 
< Fe >= (Fe5270 -t- Fe5335)/2 
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Figure 5. <Fe> is plotted against Mgb. Symbols are as Fig. 2. 
Grid lines show TMB03 models of solar abundance ratio (dashed 
lines) and 'a'— element enhancement of -1-0.3 (solid lines). Ages 
and metallicities of grid lines are as Fig. 4. The data are generally 
consistent with grid lines with 'a'-element enhancement of -t-0.3 
dex. 



galaxies (e.g. Trager et al. 2000; Proctor & Sansom 2002; 
Mehlert et al. 2003; Proctor et al. 2004b) In the following 
section we outline the procedure by which values of age and 
[Fe/H] and the 'a'-element abundance ratio are determined. 



3 DETERMINATION OF AGE AND 
METALLICITY 

For all three data sets measurements of log(age), [Fe/H] and 
'a'-abundance ratio (the derived parameters) were made 
using the x^~niinimisation technique detailed in Proctor 

& Sansom (2002) and Proctor et al. (2004a,b). The 'a'- 
abundance ratio is parameterised by [E/Fe]; the abundance 
ratio of a group of elements including C, O and Mg (see 
TMB03 for details) . Briefiy, the technique for measuring the 
derived parameters involves the simultaneous comparison of 
as many observed indices as possible to SSP models. The 
best fit is found by minimising the deviations between obser- 
vations and models in terms of the observational errors i.e. x- 
The rationale behind this approach is that, while all indices 
show some degeneracy with respect to log(age), [Fe/H] and 
[E/Fe], each index does contain some information regard- 
ing each parameter. In addition, such an approach should 
be relatively robust with respect to many problems which 
are commonly experienced in the analysis of spectral indices. 
These include poor flux calibration or sky subtraction, stray 
cosmic rays and central emission (common in LINERs and 
AGN). The method is also more robust with respect to the 
uncertainties in the modelling of both the SSPs (e.g. the 
second parameter effect in horizontal branch morphologies) 
and abundance ratio effects. It was shown in Proctor et al. 
(2004a) that, at least for Galactic globular clusters, the re- 
sults so derived are significantly more reliable than those 
based on only a few indices. 
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0.174 


-0.047 


0.445 
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0.003 









Table 1. Calibration details. For Gemini and HET data sets average errors and the difference in indices caused by the flux-calibration 
procedure (given as un-fluxed minus fluxed) are shown. Flux offsets are generally smaller than the average errors. For the SB05 data 
set average errors and the offset to the Lick system derived from observations of standard stars are shown. Lick offsets are of similar 
magnitude to the average errors. 



3.1 Reliability of the method used 

In order to further demonstrate the improvement in re- 
liability of derived parameters achieved using the 
minimisation method over the more traditional use of the 
H/3-[MgFe] index combination, we report a simple, repro- 
ducible test. First, we take the index values from a range 
of SSP models (here we use mctallicitics from -2.25 to +0.5 
dex in 0.25 dex steps and ages of 1, 1.5, 2, 3, 5 ,8, 12 and 
15 Gyr). Then we perturb each index value of each model 
SSP to represent the observational error. Here we take the 
error from a randomly generated Gaussian distribution, the 
width of which was taken from an actual SB05 observation. 
A low signal-to-noise observation was selected with Lick in- 
dex errors close to the maximum we deem suitable for age 
and metallicity determination (i.e. H/3 error >0.3 A). In this 
way a realistic model data set is generated, whose true ages 
and metallicitics arc well defined. Next, we estimate ages 
and metallicitics using the methods to be compared. 

The results of such a comparison are shown in Fig. 6. It 
is clear from Fig. 6 that the x'^-fitting of all 25 indices (with 
an rms scatter about model values of 0.07 dex or 17%) is far 
superior in accuracy to values derived from the H/3-[MgFe] 
index combination (with an rms scatter about model values 
of 0.20 dex or 58%). Furthermore, the rms scatter about 
model values when all the Balmer lines arc omitted from 
the fitting procedure (0.09 dex or 23%) is also much smaller 
than that achieved using the H/3-[MgFe] method. 

The rms in metallicity about the model values was 
0.1 dex both when all indices were included and when the 
Balmer lines omitted. For the H/9-[MgFe] method this rms 



was 0.2 dex. For [E/Fc] rms scatters wore <0.1 dex both 
when all indices were included and when the Balmer lines 
were omitted. [E/Fe] can not be estimated using the H/9- 
[MgFe] method. 

The improvement in accuracy of our x'^^minimisation 
technique over the traditional H/3-[MgFe] method could be 
understood simply in terms of the increased number of both 
age and metallicity sensitive indices employed. However, we 
have also demonstrated that excellent results are still ob- 
tained with all the Balmer lines excluded from the deriva- 
tion. This emphasises the point above that even metallicity 
sensitive indices contain some information regarding age. 
This provides one of the main motivations for our use of the 
minimisation technique. 

In the analysis reported above, we have included all 25 
Lick indices. Real observations rarely include the full wave- 
length range, and often other indices must also be omitted 
(see below). In addition, due to overlapping band definitions, 
the errors of a number of indices are correlated. We there- 
fore tested more realistic combinations of indices. These were 
based on the indices used in the analysis of the three data 
sets below, and included only 15 indices (i.e. excluding 116a, 
RSf, CNi, CNz, Ca4455, Mgi, Mgz, NaD TiOi and Ti02). 
It should be noted that the excluded indices include the ma- 
jority of those with correlated errors. The rms scatter when 
just these 15 indices were included increased to 0.1 dex (26% 
in derived age), still significantly better than those derived 
by the H/3-[MgFe] method. 
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Figure 6. Comparison of ages derived using three different com- 
binations of indices. Filled circles represent results from 
minimisation of all 25 indices. These results show a scatter of 
0.07 dex about the model values. Open circles represent results 
from x^~iiiiiiimisation of the remaining 20 indices when the five 
Balmer lines are excluded. The rms scatter about model values in 
this case is 0.09 dex. Open squares represent results from the H/3— 
[MgFc] technique. These results show an rms of 0.20 dex about 
model values. The x'^^minimisation technique is clearly superior. 
Indeed, even when all the Balmer lines arc excluded from the 
fitting procedure, the results provide a much improved accuracy 
over the H/3-[MgFe] technique. 



showed that this occurred when the index fell on a chip-gap 
interpolation (sec Fig. 1). 

- Some other individual indices were also found to be 
poor fits to the data. These were generally associated with 
sky-line residuals and vignetting (edge) effects in the data. 

The indices identified above as giving poor fits to the 
models were excluded from the fitting procedure. The av- 
erage deviations of indices to the best fit model values are 
shown in Fig. 7, where the Mgi, Mg2 and Ca4455 discrep- 
ancies are clearly visible. 

The effects of flux-calibration and the application of the 
Lick offsets of SB05 to the Gemini and HET data sets were 
also investigated. The small offsets in indices caused by flux- 
calibration (Table 1) had only small effects on the derived 
parameters. For the Gemini data a difference of +0.05±0.05 
dex was found in all three parameters. For the HET data 
there was an offset in metallicity of -^-O.liO.OS dex when 
flux-calibration was applied. Differences of -f0.05±0.05 dex 
were found in log(age) and [E/Fe]. The values of the best 
fits were found to be marginally lower for the fiux-calibrated 
data. We therefore present data for fiux-calibrated spectra 
in the following. 

Differences caused by the application of Lick offsets 
of SB05 to the Gemini and HET data sets (Table 1) also 
only caused small offsets in derived parameters (log(age): - 
0.08±0.09, [Fe/H]: 0.00±0.07 and [E/Fe]: -F0.05±0.05). The 
values of the best fits were found to be marginally lower 
when the Lick offset was not applied to the Gemini and 
HET indices. These offsets were therefore not applied to 
these data sets in the results that follow. 



3.2 Model fitting 

In the following analysis we have compared the three sets 
of observational data to two sets of SSP models; the red 
horizontal branch (RHB) models of TMB03 and the mod- 
els of TMK04. As noted previously, Fig. 4 shows that the 
TMK04 SSP models (which include a blue horizontal branch 
component) are in poor agreement with the low metallic- 
ity population in NGC 821. The models of TMB03, which 
model purely red horizontal branches, are in much better 
agreement. Given that other differences between TMB03 
and TMK04 models are small the results presented in this 
work are those derived by comparison to the TMB03 (RHB) 
models. The comparison to results derived using TMK04 
models will be summarised in Section 4. 

During the fitting procedure, a number of indices in 
each data set were found to suffer from problems that make 
them unsuitable for determining ages and metallicities (see 
Fig. 7): 

- For both Gemini and HET data sets Mgi and Mg2 
were found to be aberrant, even in the flux calibrated data 
(see Fig. 3). 

- The Ca4455 index was found to be depressed in all 
three data sets. This effect has been noted by previous 
authors in other early-type galaxies (e.g. Worthey 1992; 
Vazdekis et al. 1997; Trager et al. 1998; Proctor & Sansom 
2002). 

- In some Gemini/GMOS spectra the Fe5015 index was 
found to be a poor fit to the models. Subsequent inspection 



4 RESULTS 

The results of x'^~fittiiig of the data to the RHB models 
of TMB03 are shown in the left-hand plots of Fig. 8. The 
Gemini and HET data fitted were those that had been flux- 
calibratcd but without the application of the Lick offsets of 
SB05. Errors are estimated using 50 Monte Carlo realisa- 
tions of each data point. 

Fig. 8 shows that, generally, agreement between the 
three studies is good, despite the differing slit widths and 
orientations. Strong gradients in both ago and [Z/H] are 
evident. Fig. 8 also shows that the age gradient extends 
from the edge of the seeing affected central region out to 
log(r/rc)~0.5, with the population outside this region be- 
ing uniformly old (i.e. ^^12 Gyr). The radial [Z/H] gradi- 
ent (-0.79 dex/dex for points with -1.7<log(r/ro)<0.0) is 
extremely strong (e.g. Proctor 2002; Mehlert et al. 2003), 
suggesting that the young, central population formed from 
highly enriched gas. 

Age and metallicity determinations were also made us- 
ing the Gemini data with all the indices increased by 10%. 
This was done in order to test the impact of scattered light 
on the derived parameters in the outer regions of the galaxy 
(Section 2.1). The resultant ages for the outermost points 
were identical to those presented hero, while [Fe/H] and 
[Z/H] were increased by 0.1 dex. Given the good agreement 
between the ages and metallicities from the Gemini data and 
those from the other two studies (for which we are not able 
to estimate scattered light directly), we conclude that the 
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Figure 7. Deviations from best fit values for data compared to the models of TMB03. Mgi, Mg2 and Ca4455 can be seen to be poor fits 
in the Gemini and HET data. These indices were therefore excluded from the fitting procedure. Ca4455 is also depressed in the SB05 
data but its exclusion does not affect the fitting results. 



contamination is also <10% in all three studies. The metal- 
licity gradient given above may therefore be overestimated 
by as much as 0.1 dex/dex due to scattered light contami- 
nation. In addition, the use of TMK04 models results in a 
slight (~ 10%) reduction in the metallicity gradient (as well 
as an increase in central age of ~2 Gyr to 6 Gyr). However, 
even when both scattered light and modelling uncertainties 
are considered, the metallicity gradient remains steeper than 
-0.5 dex, leaving the qualitative results outlined below un- 
changed. 

Also indicated in Fig. 8 are estimates of the radii within 
which 5%, 10%, 20% and 50% of the galaxy's stellar mass 
resides. These radii were estimated using the luminosity pro- 
file of Ravindranath et al. (2001), an effective radius of 50 
arcsec (de Vaucouleurs et al. 1991) and the assumption of 
a constant stellar mass-to-light ratio. It is evident from Fig. 
8 that no more than ~10% of the stellar mass can have 
been involved in the star formation event responsible for 
the young central age. 

Fig. 8 also shows that the [E/Fe] values in NGC 821 are 
consistent with no gradient. However, Proctor et al. (2004b) 
expressed concerns regarding the values of [Z/H] and [E/Fe] 
that are derived directly from the SSP models. These con- 
cerns stem from consideration of the oxygen abundance ra- 
tios observed in both the solar neighbourhood, the bulge and 
the ISM of distant galaxies. 

4.1 The local abundance ratio 

Oxygen is an important element for two reasons: It is by far 
the most abundant metal - constituting ~55% of the total 
metal content of the Sun (Asplund, Grevesse & Sauval 2004). 
It is a good tracer of Type II supernovae, since its local 



Galactic abundance pattern is consistent with it having no 
other nucleo-synthetic sources (Qian & Wasserburg 2001). 
The oxygen abundance [O/Fe] in the solar immediate neigh- 
bourhood (the thin and thick discs and the halo) is approxi- 
mately -1-0.6 dex at metallicities below [Fe/H]=-1. Between 
[Fe/H]=-1.0 and [Fe/H]=0.0 [O/Fe] falls from +0.6 dex to 
0.0 dex (Chiappini, Romano & Matteucci 2003). Bensby, 
Felzing & Lundstr0m (2003) show that this downward trend 
in [0/Fc] continues at [Fc/H]>0.0, while the abundance ra- 
tios of other Q-elements (e.g. [Mg/Fe]) remain almost con- 
stant at the solar value. Fulbright, Rich & McWilUam 2004 
have recently shown a similar trend in [O/Fe] with [Fe/H] in 
Galactic bulge stars, but with [Mg/Fe] elevated with respect 
to disc stars. The low [O/Mg] ratios implied by these results 
have also been observed in the ISM of early-type galaxies 
(Humphrey & Buote 2005 and references therein). 

These observations present a challenge to current mod- 
els of Galaxy formation, and indicate that new physics is 
required (e.g. see Fulbright et al. 2004). Our concerns are 
therefore twofold: 

1) The [Z/H] values assigned to the TMB03 (or any 
other) SSP models of indices do not include oxygen. Conse- 
quently, using the up-to-date solar abundances of Asplund 
et al. (2004), we calculate that the [Z/H] values at low metal- 
licities are understated by as much as 0.2 dex, while at high 
metallicities they are overstated by as much as 0.15 dex. 

2) It is a commonly held belief that all [a/Fe] ratios 
([O/Fe], [Ne/Fe], [Mg/Fe], [Si/Fe] etc) can act as indicators 
of formation time-scale in stellar populations. Key assump- 
tions in the reasoning for this belief is that the [Fe/H] of a 
population is a good proxy for its time of formation (due 
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to Type la supernova production of Fe), while [a/H] is a 
proxy for the total amount of star formation preceding the 
formation of the population in question (due to the Type II 
supernova production of a-elements). The constant [Mg/Fe] 
and falling [O/Fe] with [Fe/H] in high metallicity disc and 
bulges stars of our Galaxy clearly indicates a breakdown in 
one or more of these assumptions. We must therefore either 
abandon the idea that the [Fe/H] of a population is a proxy 
for its time of formation (and the belief that [a/Fe] ratios 
can act as indicators of formation time-scales), or we must 
accept that not all [a/Fe] ratios can act as indicators of 
formation time-scale in stellar populations. 

As a result of these considerations we have applied cor- 
rections to the values of [Z/H] to include the effects of the 
oxygen abundance variations in the stars used to construct 
the SSP models. Both corrected and uncorrected values are 
presented in the results that follow. 

The correction to [Z/H] was carried out by applying an 



offset (AZ) to the values derived directly from the TMB03 
models using the up-to-date solar abundances of Asplund 
et al. (2004), the [O/Fe] trend outlined previously. AZ was 
given by: 

AZ=-h0.175 (for [Fe/H]<-1.0) 

AZ=-0.175[Fe/H] (for -1.0<[Fe/H]<0.0) 

AZ=-0.33[Fe/H] (for [Fe/H]>0.0) 

Rather than correct [E/Fe] for local trends in [O/Fe], 

wc have chosen to represent the 'a'-element abundance ra- 
tio by [O/Fej. Both [E/Fe] and [O/Fe] are presented in the 
results that follow. [O/Fe] was calculated using: 

[O/Fe]=[E/Fe]-F0.35 (for [Fe/H]<-1.0) 
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Figure 9. Colours predicted by age and metallicity from TMB03 
SSP models and the colours of TMK04. Symbols arc as Fig. 2. 
The radial trends in colours from Goudfrooij et al. (1994) are 
also shown (lines). However, offsets (given top-right in each plot) 
have been applied to the Goudfrooij et al. data in order to match 
our results. These are consistent with the effects of extinction 
but other contributory factors may be present, such as modelling 
uncertainties in the SSP models. 



[O/Fe] = [E/Fe]-0.35[Fo/H] (for -1.0<[Fe/H]<0.0) 

[O/Fe]=[E/Fe]-0.60[Fe/H] (for [Fe/H]>0.0) 

It is important to note that the expressions given above 
assume that the [O/Fe] with [Fc/H] trend in the stellar pop- 
ulation under consideration is the same as that in the solar- 
neighbourhood. This assumption is supported by the fact 
that in [O/Fe] with [Fe/H] trends of stars in the bulge of 
our own Galaxy show an identical trend to the stars in the 
disc (Bensby et al. 2003; Fulbright et al. 2004), even though 
the [Mg/Fe] of the two populations differ significantly. We 
shall refer to the correction to [Z/H] and the transform from 
[E/Fe] to [O/Fe] as the 'local abundance ratio correction' 
(LARC). 



licity, colours were estimated from the models of TMK04 
using the age and [Z/H] values shown in Fig. 8. Colours 
for age and [Z/H] values derived directly from the TMB03 
SSP models are shown in the left-hand plots of Fig. 9, while 
colours from LARC values are shown in the right-hand plots. 
These estimates are compared to the photometric observa- 
tions of Goudfrooij et al. (1994). It was necessary to apply 
significant offsets (given in the plots) to the Goudfrooij et 
al. (1994) data in order to obtain a match between our es- 
timates and the observations. The sense of these offsets are 
consistent with the effects of internal extinction, although 
other systematic effects (e.g. SSP modelling errors) may be 
present. However, once this offset is applied the match be- 
tween observations and data are extremely good, i.e. we re- 
cover the relatively shallow colour gradients of Goudfrooij et 
al. (1994) extremely well, despite the strong age and metal- 
licity gradients present in our spectroscopic results. 



5 DISCUSSION 

The central luminosity-weighted age of ~4 Gyr found in this 
work is consistent with the young ages found in some pre- 
vious studies. However, as noted in the introduction, other 
studies have found old ages. We find a natural explanation 
for these results when we consider the strong age gradient 
and the varying aperture size between studies. The results 
presented here for the central regions are for small regions 
(0.5x0.5 arcsec, 1.0x0.3 arcsec and 2.1x1.1 arcsec for the 
Gemini, HET and SB05 data sets respectively). Denicolo et 
al. (2004), who obtain a similarly young age of 4.2 Gyr, used 
an aperture 1.5x12.0 arcsec, while the Trager et al. (2000) 
age of 7.7 Gyr was for a 2.0x12.0 arcsec. The two studies 
that found old ages also used large apertures - Vazdekis et 
al. (2004) quote an age of 12.5 Gyr within an aperture of 

I. 6x10 arcsec, while Caldwell et al. (2003) find an age of 

II. 5 Gyr using a 3x5 arcsec slit. A plot of effective aper- 
ture size (using the j0rgensen, Franx & Kjaargaard 1995 
definition) the against the measured age shows a correla- 
tion, albeit with significant scatter. We therefore conclude 
that the presence of the age gradient and the varying aper- 
ture shapes/sizes between the various studies plays a part 
in causing the range of ages in the literature. 



There are two significant effects of the local abundance 

ratio correction on our results for NGC 821. The effect on 
the [Z/H] gradient is to reduce it somewhat (to ~-0.60 
dex/dex), although this is still an extremely strong metal- 
licity gradient. The effect on the slope of the enhancement 
parameter (now [O/Fe]) is quite profound, with a gradient 
evident in [O/Fe] that was not present in [Mg/Fe]. The gra- 
dient is such that the young, central population is signif- 
icantly less enhanced in oxygen with respect to iron than 
the older, outer regions. This is, at least qualitatively, more 
in line with expectations, given that we might expect the 
young, metal-rich population in the centre to have formed 
from gas that had been enriched by Type la supernovae over 
a longer time-scale than the old population in the outer re- 
gions of the galaxy. 

As a further test of the derived values of age and metal- 



We next turn our attention to possible star formation 

histories in NGC 821. The young central age clearly sug- 
gests relatively recent star formation, although there is no 
evidence for current star formation. It is therefore appropri- 
ate to consider the source of the gas that fuelled this star 
formation. By making some simple Eissumptions regarding 
the relative luminosities of young and old populations we 
estimate that the mass involved in the more recent burst 
to be between a few times ~10*Mq (assuming the burst to 
be 1 Gyr old) and a few times lO^^M© (assuming the burst 
to be 4 Gyr old). We therefore constrain the mass involved 
in the recent starburst to be bctwcou ^0.2% and ~10% of 
the overall galactic mass. These values above are in stark 
contrast to the current mass of dust (~2.5x lO'^M©; Bettoni 
et al. 2003) and HI gas (<10**Mq; Georgakakis et al. 2001). 
It is therefore clear that the gas involved in the recent burst 
was either almost completely consumed in star formation or 
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was ejected from the galaxy as a wind. 

Next, wc consider a number of sources for the gas that 
fuelled the recent starburst; a major merger, a minor gas- 
rich merger and in-situ gas. 

The first consideration is the extremely strong metal- 
licity gradient present in NGC 821. Such a strong gradient 
is inconsistent with merger models which predict small gra- 
dients due to the 'washing-out' of the ambient gradients in 
the turbulent mixing of the merger (Kobayashi 2004) . The 
gradient in NGC 821 (~ 0.7 dcx/dcx) is therefore signif- 
icantly greater than is achieved in merger models. In ad- 
dition, NGC 821 lacks the tidal tails, shells and plumes 
expected to persist for ~3 Gyr after a major merger (e.g. 
NGC 1700; Brown et al. 2000). Thus, a recent major merger 
seems unlikely. 

An alternative explanation is that of the accretion of a 
small, gaseous, satellite galaxy. However, in this scenario, the 
extremely high central metallicities are hard to achieve since 
the small accreted galaxy would, presumably, have contained 
low metallicity gas. 

The final scenario we consider is that of a central burst 
of star-formation fuelled by in-situ gas from the galaxy itself. 
This could have come from the smooth accretion of cooled 
halo gas onto a disc as predicted by hierarchical models (e.g. 
Meza ct al. 2003). It is therefore possible that the recent 
starburst was associated with gas from the embedded stel- 
lar disc (Goudfrooij et al. 1994; Scorza & Bender 1995; de 
Souza et al. 2004). However, Scorza & Bender (1995) show 
the disc to be highly inclined, with an inclination of ~72°. 
Therefore, since the young ages are seen to extend to the 
same relatively large radii on both major and minor axes 
(Fig. 8), it is unlikely that the recent starburst was confined 
to the disc. 

Regardless of the source of the gas within the galaxy, 
we would expect the younger population to have 'a'-element 
abundances depressed with respect to the outer regions. In 
other words, we would expect the time-scale of star forma- 
tion to be longer in the young, high-metallicity, central re- 
gion of NGC 821 than in the old, low-metallicity outer re- 
gions. This is clearly not the case for the [E/Fe] values de- 
rived directly from TMB03 models. It is the case, however, 
for our estimates of [0/Fe]. We are therefore able to con- 
struct a recent star formation history for NGC 821 that is 
consistent with both spcctroscopically determined ages and 
metallicities and the photometric data by assuming that be- 
tween 0.2% and 10% of the galaxy's mass was formed 1 Gyr 
to 4 Gyr ago, respectively, from the galaxy's own gas. The 
paucity of gas remaining within the galaxy favours the lower 
end of this mass range and younger ages for the recent star- 
burst. This, then, represents our favoured scenario. 



6 CONCLUSIONS 

We have presented the analysis of Lick indices from three 
separate long-slit spectroscopic observations of the nearby 
elliptical galajcy NGC 821. The three data sets have been 
shown to present a consistent picture of the stellar popu- 
lation within one effective radius. Strong gradients in both 
age and metallicity were observed. The old (> 10 Gyr), low 
metallicity ([Fe/H]~-1.0) population in the outer regions 



of the galaxy (r^^ Itc) exhibit no sign of the blue horizon- 
tal branch morphologies present in some Galactic globular 
clusters of the same age and metallicity. We have shown that 
the large range of ages in the literature is probably the re- 
sult of the range of aperture sizes and the strong central age 
gradient. 

We therefore find the central region of NGC 821 to have 
undergone a burst of star formation between 1 Gyr and 
4 Gyr ago. This conclusion is based on the 4 Gyr luminosity- 
weighted age found for the integrated spectra of the central 
region and the observation that the galaxy, which shows no 
current star formation, has also passed through the post- 
starburst phase (which persists for approx. 1 Gyr) . The old 
ages in the outer regions indicate that the burst was cen- 
trally concentrated and we have shown that it can not have 
involved more than 10% of the galaxy's stellar mass. 

We find the most likely scenario to be a very recent 
(~1 Gyr), low mass (~0.2% of total galaxy mass) starburst 
in the in-situ gas, perhaps triggered by the accretion of a 
small (gas poor) satellite galaxy or a tidal interaction. In 
this scenario, the ambient gas in NGC 821 was consumed in 
a centrally concentrated burst of star formation. The burst 
would have had a high [Fe/H], but low [O/Fe], due to the 
long (>10 Gyr) time-scale of star formation in the galaxy. 
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